Quantum Error Correction with Uniformly Mixed State Ancillae 
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It is often assumed that the ancilla qubits required for encoding a qubit in quantum error correction 
(QEC) have to be in pure states, 1 00 . . . 0) for example. In this letter, we introduce an encoding 
scheme avoiding fully correlated errors, in which the ancillae may be in a uniformly mixed state. We 
demonstrate our scheme experimentally by making use of a three-qubit NMR quantum computer. 
Moreover, the encoded state has an interesting nature in terms of Quantum Discord, or purely 
quantum correlations between the data-qubit and the ancillae. 



A quantum computer is vulnerable against environ- 
mental noise and it must be protected by one way or 
another. Quantum error correction (QEC) is one of the 
most successful approaches to this end [l|. Despite this 
great success, QEC requires expensive resources, or an- 
cillae that are usually assumed to be in pure states 0, Q . 
However, it is not yet proved that ancillae in uniformly 
mixed states are useless. We extend previous works [4] 
and show an encoding scheme robust against fully corre- 
lated noise in which all the ancillae can be in uniformly 
mixed states. The encoded state has an interesting na- 
ture in terms of Quantum Discord [1| , or purely quantum 
correlations between the data-qubit and the ancillae. Our 
QEC scheme also provides an example of Deterministic 
Quantum Computation with 1-Qubit (DQC-1) [1, 0j- 

Suppose we have a single qubit in an arbitrary state 
pi, which we want to protect from noise. We introduce 
some additional qubits (ancillae) in order to protect the 
first qubit and suppose that all the qubits suffer from the 
same noise. Such a noise is called fully correlated and it 
may happen when the dimensions of the quantum com- 
puter are microscopic compared with the wavelength of 
external disturbances. Noiseless subsystem (NS) [a-tll| 
and decoherence free subsystem (DFS) (l2rll5j are well 
known strategies to protect a system from such fully cor- 
related noises [THUS]- These schemes, however, require 
ancillae in pure states and thus they are expensive. 

In the following, we show that it is indeed possible to 
devise a cheaper QEC scheme employing ancillae in the 
uniformly mixed state. Let 



Pi 



a /2 + (n x ,n y ,n z ).(a x ,a y ,a z )/2 (1) 



be the state of the qubit to be protected. Here Co is 
a unit matrix of dimension 2, n = {n Xl n y ,n z ) is the 
Bloch vector, and <ii is the ith component of the Pauli 
matrices. We introduce two ancillae in uniformly mixed 
states, whose Bloch vectors are 0. The initial state of the 
three-qubit system is thus a tensor product state p\ ® 
(cto/2)® 2 . The unitary encoding operator ?7e transforms 
the tensor product state to an entangled state p 3 . If the 
state of the system is again p^ even after the action of 
noises, a unitary recovery operator, {J R = [/jj,, transforms 
p~3 back to the initial tensor product state pi (g> (cr /2)® 2 
and pi can be recovered after tracing over the ancilla 



states. 

It is highly counterintuitive that a QEC scheme works 
with ancillae in uniformly mixed states. The trick is that 
the uniformly mixed state (cro/2)® 2 is rewritten as 



j {\n 2 ,n' 2 ){n 2> n' 2 



- n 2 ,-n 2 )(-n 2 ,-n 2 \ 
-n 2 ,n' 2 ){-n 2l n' 2 \ + |n 2 , -n' 2 ){n 2 , -n' 2 \ ), (2) 



where n 2 and n' 2 are arbitrary Bloch vectors ( | n.2 1 = 
\n' 2 \ = 1) and \n 2 ) and \n' 2 ) are pure states correspond- 
ing to n 2 and n' 2 , respectively. If a QEC scheme works 
with arbitrary pure ancilla states, the superposition prin- 
ciple of quantum mechanics guarantees that ancillae in a 
uniformly mixed state do work as well. 

A more formal description is given as follows. Suppose 
we have a single qubit in a state pi, which we want to 
protect from noise operators {a x ,a y , a z }. To this end, we 
introduce two additional qubits, which may be in an ar- 
bitrary state p 2l and apply a suitable encoding operator 
Ue onpi®p 2 to obtain a codeword p~3 = Ue {pi ®P2)U E - 
We introduce the fully correlated error channel $ repre- 
sented by 
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®3 



(3) 



Here 



where X = af 3 .X 1 = af 3 ,A 2 = <rf 3 ,X 3 = a 
Pi > is the probability with which an error operator Xi 
acts on p3 and we assume X}j=o Pi 
an encoding operator Ue satisfying 

U E XiU E - cr Cg) Mi 



1. Suppose there is 



(4) 

for i — 1,2,3. Then, Ue defines the QEC scheme that 
we are seeking. We can show that 



U R ${U E {pi®P2)U E )ul = 



i=0 

Pi®p' 2 , 



(5) 



where p' 2 = J^iPiMiP^Mi ■ This proves that, after de- 
coding, the error channel $ affects only p 2 but not p\. 

There are infinitely many choices of Ue but care- 
ful inspection of the error operators reveals that Ue = 



was shown for this scheme that 




U E U r = U, 



FIG. 1: Encoding circuit Ue, error channel <£> and recovery 
circuit Ur in the simplest case. 



C^CNOT3i £^cnoti2 is the simplest choice [!, fl8j |. The 
4x4 matrices {Mi} are obtained by direct calculation as 
M = af 2 , Mi = erf 2 , M 2 = -o y ®a x and M 3 = a z ®a . 
Figure Q] shows the encoding circuit Ue, the error channel 
$ and the recovery circuit Ur. 

Let tSi and S 2 be the state spaces of the data qubit and 
the ancillae, respectively. The set of the encoded states 
is a subset of the total state space £3 
of the three-qubit system. How our QEC scheme works 
is summarized in Fig. [5] 

The extreme case of the uniformly mixed state p 2 = 
erf 2 / 4: is worth analyzing separately. When the data 
qubit is in a pure state, this provides an interesting ex- 
ample of DQC-1 Moreover, our QEC scheme equally 
works for a data qubit in a mixed initial state. It is read- 
ily found that 

/ / n-® 2 \ \ rr® 2 

U a * \U B (p! ® Ulj U R =pi®\. (6) 

This shows that pi eg) cr® 2 /4 is a fixed point of this oper- 
ation for any p\ G <Si. 

Let us compare our QEC scheme with the NS encoding 
scheme discussed in This NS encoding scheme em- 
ploys three qubits to encode a logical qubit robust against 
any noise of the form W® 3 , where W is an arbitrary el- 
ement of the 2-dimensional representation of SU(2). It 
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FIG. 2: Encoding operator Ue maps a tensor product state 
pi <S> P2 € Si ® 52 to a codeword p3 € S3. The error channel 
$ maps pz to p' 3 within the code space 53. As a result, the 
recovery operator Ur maps p' 3 to pi® p2, restoring the data 
qubit. 



u' R *' (u' e (pi ® |0> <0| ® p a )t/^) c/^ 1 

3 

= pi® \0)(0\ Q^ptUiPaU}, 
i=0 

where L/^ and U' R are the encoding and the recovery op- 
erators, respectively. {Ui} = {a , e la<Ta! , e 2 ' 3 ' 7 " , e I7CT " } is 
the set of error operators and <&' is defined by an expres- 
sion analogous to Eq. ([3]). For this scheme, the initial 
ancilla state is required to be of the form |0)(0| ® p a . In 
contrast, although the error operators avoidable with our 
scheme are restricted within a subset of those avoidable 
in |4(, our proposal has the remarkable advantage that 
any initial ancilla state P2 can be employed for successful 
QEC. 

We will discuss quantum discord (hereafter, abbrevi- 
ated as QD) introduced in Q in order to analyze an- 
other aspect of our scheme. QD is a measure of non- 
classical correlations between two subsystems of a quan- 
tum system. Surprisingly enough, it was found that QD 
may be non-vanishing even in the absence of entangle- 
ment and that in fact, there are useful quantum algo- 
rithms that work with little or no entanglement called 
DQC-1 ism. In other words, separability alone does 
not imply the absence of a quantum nature of the state. 

The left and the right QDs of our encoded state 

1 / 

P3 = o (. fT ° ® a o ® cr + n x a x ®a x ®a 
o 

+ n y a y ® o x ® a z + n z a z ® a ® a z ) , (7) 

are defined respectively as [!, [2(| 

2?(2:1) = S(p 1 )-S(p 3 ) + S(2\l), 
P(l:2) = S(p 2 )-5(p 3 )+S(l|2). 

Note that they arc not necessarily equal to each other Q . 
Here, S(p) = — Tr(plog 2 p) is the von Neumann entropy 
of a density matrix p. The density matrices pi and p 2 
are obtained by tracing over the ancillae and data-qubit 
states, respectively. We define a projective measurement 
by a complete set of two orthonormal vectors {| ± m)}, 

which define | ± m)(±m| = — - — with \m\ = 1. 

Let us define the conditional entropy by Sn± m \\ (2|1) = 

y^ J P±S(P2\±), where p 2 \± = (±m\p 3 \ ± m)/p± and 
± 

p ± = Tr(±m|p 3 | ± m). Then, S(2\l) is defined as 

min 5/i± m )}(2|l). 5(1|2) is defined similarly. 
|m|=l 

The explicit form of the ancilla state after the mea- 
surement of the data qubit with a basis {| ± m)} is 

1 . 

P2\± = I WO ® Co ± n x m x a x ® cr 

± n y m y a x ® a z ± n z m z <To ® <j z ) . 
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FIG. 3: Example of T>{\± m )}(2 : 1) when the initial state of 
the data qubit has the Bloch vector n — x. (r, 9 
(2?{ | ± m )} (2 : l),S m ,0m) is plotted, where m and <f> m specify 



m = (sinf 



0- ^{|±™>}(2 : 1) = 



for m — ±x. Therefore, 23(2 : 1) vanishes for n — x. 

The corresponding conditional entropy is 
1 8 

£{|±m>}(2|l) = 2--2^(l + n-nj-)log 2 (l + n-nj), 
i=i 

where = (±m x , ±m y , ±m z ) are all eight combinations 
of three ±. 

As an example, we show T>^ m ^(2 



1) = S(/5i) - 

S(p3) + S{|± m )}(2|l) for the initial state p\ with 
n = i as a function of m = 
m sm4> m ,cos9 m ) in Fig. El Here, x,y 



the Bloch vector 
(sin 9 m cos (j) m , sin ( 
and z are the unit vectors along the x-, y- and z-axes, 
respectively. Note that when m = x, Pru^i = 0. 
Therefore, 23(2 : 1) = 0. The quantum discord 23(2 : 1) 
as a function of n = (sin 9 cos 4>, sin 9 sin (f>, cos 9) is shown 
in Fig. SJ 




FIG. 4: Quantum discord 23(2 : 1) as a function of the ini- 
tial state of the data qubit parameterized by the Bloch vector 
n = (sin 9 cos <f>, sin 9 sin <j>, cos 9). See Eq. {TJ . Coordinates 
(r = 25(2 : 1), 9, <f>) depict QD as a function of 9, <j>. 23(2 : 1) 
vanishes when n = (±1, 0, 0), (0, ±1, 0), (0, 0, ±1). The func- 
tion 23 (2 : 1) takes the maximum value (3/4) log 2 3— 1/2 when 
n = (±1, ±1, ±1)/V3. 



Although extensive optimization is necessary to eval- 
uate QD in general, some initial states satisfying 
23 (2 : 1) = are easily obtained by carefully in- 
specting the structure of /5 3 . Let us consider the case 
n = x, for example. In this case, p^ is reduced to 
| (<7q <g) <to <E) <7o + <t x ® &x ® &o) , which contains only <tq 
and o x for the data qubit. Therefore, it is reasonable to 
employ \x) as a candidate for \m) for obtaining S'(2|l). 
With this choice, is reduced to a block-diagonal form 

p 3 (n = x) = ^2\±x){±x\ <g> (p±p 2 \±) 



and 23(2 : 1) = is readily obtained [5H20j. 

23(1 : 2) = for an arbitrary initial state can be proved 
similarly. p 3 contains only uq ® o"o, <J X ® <jq, <jo®<j z , and 
c x 65 &z for the ancilla qubits. Therefore, we take 



|n±±)(n 



± X, ±z)(±x, ±z\ = ® 



as a complete set of four unit vectors that determine the 
projective measurement on the ancillae, although there 
are many other possibilities. p 3 is rewritten as 



P3 



£(P±±Pl|±±) 
±± 



|n±±)(n± ± |. 



When the data qubit and ancillae are rearranged, the 
density matrix is rewritten as a block-diagonal form and 
thus 23(1 : 2) = is immediately obtained. 

According to the classification introduced in (2ll. |22T|. 
vanishing 23(1 : 2) implies that our encoded state has 
a quantum-classical correlation. Furthermore, in case 
23(2 : 1) also vanishes, p% has a product eigenbasis as we 
have shown above, and the encoded state has a classical- 
classical correlation, or, in other words, is (properly) clas- 
sically correlated. 

When the ancillae are pure, we find 23(2 : 1) = 23(1 : 
2), which is nothing but the entanglement entropy. For 
example, 

23(2 : 1) = 23(1 : 2) 

= 2 - (1 - n z )log 2 (l - n z ) - (1 + n 2 )log 2 (l + n z ), 

when p 2 = \z){z\ <g> \z)(z\. 

We demonstrate our QEC scheme with a NMR quan- 
tum computer. We employ a JEOL ECA-500 NMR spec- 
trometer [23[, whose hydrogen Larmor frequency is ap- 
proximately 500 MHz. We employ a linearly aligned 
three-spin molecule, 13 C-labeled L-alanine (98% purity, 
Cambridge Isotope) solved in D 2 0. 

We simplify the quantum circuit shown in Fig. [T] by 
taking into account the fact that the phases of states 
are not independently observed in a NMR quantum com- 
puter. Both the encoding and the decoding require only 5 
pulses including refocusing pulses, taking approximately 
25 ms. 



y 1 - a y/r- b y 1 c 




Error Operator Xo Xi V2 
Fe((T 0) X) 0.68 0.73 0.75 
Tr(X) 1.03 1.00 1.00 

FIG. 5: Visualization of error correction performances. The 
surface of the Bloch sphere is mapped onto the surfaces in (a), 
(b) and (c) corresponding to three different error operators, 
Xq,X\,X2, respectively. See the text for details. Entangle- 
ment fidelities F e (ao, M) and traces Tr(.M) are summarized 
in the table. A4 represents a map which is determined by the 
encoding, error, and recovery processes depicted in Fig. [T] 

The density matrix of the thermal state is well approx- 
imated by 

P = (<V2)® 3 

e / \ 

o 

where e ~ 1CT 6 . Since (cr /2)® 3 is not visible in NMR, 
the density matrix of the thermal state is considered as 
a pseudo-pure state for DQC-1. 

We perform three sets of experiments, in which we set 

( a ) {Pi} = (1)0)0, 0) : no error 

(b) { Pi } = (0,1,0,0) :Xx error 

(c) {k} = (0, 0, 1, 0) :X 2 error, 

respectively, in Eq. We do not need to examine the 
X3 error separately since X3 = 1X2 X\. Each set starts 
with 4 different initial states in order to apply quantum 
process tomography 24] . The results are summarized in 
Fig. [SJ Although the surfaces are distorted, it is clear 



that our QEC scheme indeed eliminates the effects of the 
fully correlated noises. The table in Fig. [5] summarizes 
the entanglement fidelities. 

It is noteworthy that one-qubit gate operations on 
the logical qubit take simple forms. Let V be a one- 
qubit gate acting on the logical qubit. Then its ac- 
tion on the physical qubits is obtained by simplifying 
Ue{V ®<Jo®&o)Ur. For the simple gates V = a Xl o y and 
a Zl the corresponding operations on the physical qubits 
are a x ® a x ® a® , a y ® a x ® a z , and a z ® cto ® cr z , re- 
spectively. Note that these operators satisfy the ordi- 
nary 5u(2) algebra. It is easy to obtain more general 
gate operations acting on the physical qubits by sim- 
ply exponentiating these operators, e.g., e -i- a < J ^®< J ^®' J o 
implements V = e' 10 " 7 - . Note that e -^°y®°*®°z = 

e -iTr((T z ®(To'2icr i )/4 e -i / 9cr x (g)<T :c (gicro e i7r(<T z (8(To®(T z )/4 j„ fag case 

of V = e _J/3<T!/ . From these operators, we can understand 
how the information of the data qubit is distributed in 
the encoded state. We note that direct operations on 
logical qubits in DFS/NS were discussed in [jf| . 

In summary, we demonstrated a quantum error correc- 
tion scheme avoiding fully correlated errors, in which the 
ancillae can be in a uniformly mixed state. Our results 
pave the way to new applications of DQC-1 to quantum 
computing. The analysis of quantum discord reveals that 
our encoding creates an interesting quantum correlation 
between the data qubit and the ancilla qubits; our en- 
coded state has a quantum-classical correlation in gen- 
eral and has a classical-classical correlation when both 
left and right quantum discords vanish. We anticipate 
further progress both in the understanding of quantum 
correlations and the development of QEC schemes. Our 
QEC scheme admits simple one-qubit gate operations on 
the encoded qubits. 
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